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Membrane budding and wrapping of particles, such as viruses and nano-particles, play a key 
role in intracellular transport and have been studied for a variety of biological and soft matter 
systems. We study nano-particle wrapping by numerical minimization of bending, surface tension, 
and adhesion energies. We calculate deformation and adhesion energies as a function of membrane 
elastic parameters and adhesion strength to obtain wrapping diagrams. We predict unwrapped, 
partially-wrapped, and completely-wrapped states for prolate and oblate ellipsoids for various aspect 
ratios and particle sizes. In contrast to spherical particles, where partially-wrapped states exist 
only for finite surface tensions, partially-wrapped states for ellipsoids occur already for tensionless 
membranes. In addition, the partially-wrapped states are long-lived, because of an increased energy 
cost for wrapping of the highly-curved tips. Our results suggest a lower uptake rate of ellipsoidal 
particles by cells and thereby a higher virulence of tubular viruses compared with icosahedral viruses, 
as well as co-operative budding of ellipsoidal particles on membranes. 



I. INTRODUCTION 

Budding and subsequent vesiculation of lipid bilayer 
membranes is essential for transport in biological cells 
[TH2] . Biomembranes usually consist of a variety of lipids 
and proteins, therefore buds can be induced by lipid do- 
main formation [1HZ], by membrane spontaneous curva- 
ture [SHU], or by a combination of both [TB-15J. For 
particle wrapping, in addition to the membrane proper- 
ties, the particle shape and the adhesion strength have 
to be taken into account [HH2D] . Biological examples are 
viral budding [2U [22] and — for designing efficient drug 
delivery systems and other nano-engineered techniques 
for medical diagnosis — the uptake of small particles 
by cells [T71 HOI 1231 Hi]- Elongated viruses have been 
found to form patterns on cell membranes |25j , and ellip- 
soidal nano-particles are used for drug delivery [26] and 
as markers [27J [2E] in cell biology. 

Techniques to fabricate nano-particles of different 
shapes and sizes are well established, thus wrapping has 
been studied experimentally using both vesicles and cells. 
To gain insight into the basic mechanism of cellular up- 
take, the role of shape and size has been investigated 
both in vitro QI1 El [201 123] and as well as in vivo [H- 
|3"T] studies. 

Describing the membranes as a mathematical surface 
with appropriate curvature-elastic constants, continuum 
models predict wrapping as function of bending rigid- 
ity k, spontaneous curvature cq, surface tension a, and 
for lipid domain formation also the line tension 7 at the 
domain boundary. Bending rigidity and surface tension 
oppose wrapping, whereas the adhesion strength w fa- 
vors wrapping. For example, for a completely wrapped 
sphere of radius R and for a lipid bilayer without surface 
tension, an adhesion energy gain of — 4irR 2 w is opposed 
by a bending energy cost of 8ttk. Complete wrapping oc- 
curs if the adhesion strength exceeds w* — 2k/ R 2 , while 
for smaller values the sphere remains unwrapped. 

Wrapping of spherical particles has been studied sys- 
tematically using continuum membranes and is well un- 



derstood. A particle wrapped by an infinitely large pla- 
nar membrane without surface tension is fully described 
by the simple calculation above. Neglecting a surface 
or volume constraint, also the energy of a vesicle that 
wraps a particle has been obtained analytically [32] . For 
membranes with surface tension, the deformation energy 
of the membrane can be calculated using approximate 
models [23 [33] or shape equations that are evaluated nu- 
merically [34 36 . Similarly, the deformation energy for 
wrapping of infinitely long cylinders has been calculated 
[571 155] . which is qualitatively different from wrapping 
of spheres. Regarding the example given above, the free 
membrane around the sphere forms a catenoid without 
any bending energy cost, whereas for a cylinder wrapped 
by an infinite planar membrane, the deformation energy 
of the free membrane has to be taken into account. 

In this article, we investigate wrapping of ellipsoidal 
particles by homogeneous membranes without sponta- 
neous curvature. In section [ill A[ we calculate the wrap- 
ping energy as function of the wrapping fraction for 
membranes with and without surface tension. We ob- 
tain phase diagrams that show a non-wrapped state, 
a partially-wrapped state, and a fully-wrapped state, 
see section IIIIBI While the transition between the un- 
wrapped and the partially wrapped state is continuous, 
the partially-wrapped state is separated from the fully- 
wrapped state by an energy barrier. In section [III C[ we 
characterize the energy barrier and the hysteresis that is 
found for the transition between the partially-wrapped 
and the fully- wrapped state. Finally, in section [III D| we 
discuss the role of the shape for the wrapping energy. 

II. THEORETICAL MODEL AND METHODS 

A. Continuum membrane model 

Using the continuum membrane model, our system is 
constructed in order to include the minimal ingredients 
required to characterize wrapping of ellipsoidal particles. 



FIG. 1. Deformation profile of membrane adhering to a rigid 
nano-particle. The free membrane around the particle is la- 
beled a and the membrane adhered to the particle /3. 

The uptake process can be understood as a competitive 
tug of war between the elastic deformation energy and 
the contact interaction between particle and membrane. 
The elastic deformation energy £^ e f of the lipid bilayer 

is US iD] 

S dei = [ dS [2k(H - c ) 2 + RK + a] , (1) 
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which is an integral over the entire membrane surface 
area, A to t- The membrane shape is described by the 
mean curvature, H = (ci + C2)/2, and the Gaussian cur- 
vature, K = c\C2] c\ and c-i are the principal curvatures 
of the membrane. The surface tension a is the Lagrange 
multiplier conjugate to the excess area that is generated 
due to wrapping on the particle relative to the unwrapped 
flat membrane. Because we assume a symmetric lipid bi- 
layer, we use Co = 0. The integral over the Gaussian 
curvature with the constant saddle splay modulus re is 
determined by the topology of the membrane (and by 
the geodesic curvature at a boundary). In our case of 
an infinite planar membrane, the integral is a constant 
during the wrapping process. Thus, the total energy for 
a membrane-particle wrapping complex is 

£tot = / dS [2kH 2 +a]-w f dS . (2) 

The contact interaction with adhesion strength w is pro- 
portional to the membrane area ^4 a d adhered to the par- 
ticle, see Fig. [I] 

Our continuum model in Eq. ^ is applicable for par- 
ticle sizes that are larger than a few times the thickness 
of a lipid bilayer, which is about 5 nm. For particle sizes 
smaller than \J a both bending energy and surface ten- 
sion contribute, for larger particle sizes surface tension is 
dominant [35] . A characteristic length scale for each sys- 
tem is the particle size a, typically 20 — 100 nm and a 
characteristic energy scale is the bending rigidity re, typ- 
ically 10 — 100 kftT. We therefore describe our system in 
terms of dimensionless parameters, which we choose to 
be consistent with those in Rcfs. 1551 and 131)1 This gives 

(3) 

with the reduced energy E = £ to t/(irK), the reduced 
surface tension a = era 2 /re, and the reduced adhesion 
strength w = 2wa 2 / 're. 




FIG. 2. Membrane deformation for membrane-particle sys- 
tems that have been discretized using triangulated surfaces. 



B. Energy minimization 

There are three approaches to calculate the membrane 
deformation for particle wrapping |41j . (i) Solving the 
Euler-Lagrange equation: the Euler-Lagrange equation 
is obtained from Eq. ^ using variational calculus; it 
is a fourth-order partial non-linear differential equation 
[32] and a general solution does not exist. Exploiting 
symmetry and choosing an appropriate parametrization 
[4TJI |4*3l [33] , such as the arc-length parametrization pro- 
posed in Ref. [331 one obtains a set of second order non- 
linear ordinary differential equations. This method has 
been employed both for spherical particles in Refs. l34l — 
[31)1 and H5] and for infinite cylinders in Refs . 1381 and Hoi In 
the weak-wrapping limit, the Monge parametrization and 
the approximation for almost planar membranes can be 
used |34H37| . (ii) Exploiting a variational approach: min- 
imum energy shapes are found using appropriate func- 
tional parametrization, e. g. spherical harmonics [46] . 
Cassini ovals [39], and Fourier functions [47] • The set of 
parameters obtained from energy minimization describes 
the membrane deformation. This method has been used 
for wrapping in Ref. 48 and can treat non-axisymmetric 
shapes as well, (iii) Using triangulated membranes: the 
method offers a high flexibility and has been used to 
study both minimal energy shapes 49, 50] and systems at 
finite temperature |5 1H54| . The membrane is constructed 
by small triangles and discretized expressions for the de- 
formation and interaction energies are used. Triangu- 
lated membranes provide a huge amount of freedom re- 
garding the choice of shape and local properties of the 
membrane and has been used to study wrapping of single 
as well as multiple particles [54 - 56] . We employ "Surface 
Evolver" [3H] for our calculations, a finite-element code 
for surfaces formed of vertices, edges and facets; the dis- 
cretization of the bending energy is discussed in appendix 

m 



C. Wrapping energy calculations 

Deformation energies are calculated such that the 
membrane wraps the particle starting from the lowest 
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curvature region progressively. Therefore, ellipsoidal par- 
ticles are oriented with their major axis aligned parallel 
to a membrane patch that is enclosed by a circular wire 
frame of radius 20 a, see Fig. [2j This orientation is ener- 
getically favourable until about half wrapping compared 
with the perpendicular orientation. It ensures a max- 
imum gain of adhesion energy and a minimal cost due 
to bending and surface tension contributions. We char- 
acterize the ellipsoids by their polar radius a and their 
aspect ratio b/a. Prolate ellipsoids correspond to b/a > 1 
and oblate ellipsoids correspond to b/a < 1. 

From numerical calculations for ellipsoids with parallel 
orientation, we obtain the deformation energy as function 
of the fraction of the particle area that is wrapped, ^4 a d/A 
[57j . Snaphots of partially- wrapped particles with differ- 
ent wrapping fractions are shown in Figs. [3] (b), [4] (b), 
and |4|(d) 58J. The energy of the membrane relative to 
the unwrapped state is denoted by AE, which gives the 
wrapping energy cost AE = AE/ttk in reduced units. 
To calculate small energy differences and derivatives with 
sufficient accuracy, this deformation energy profile is fit 
to the generalized logistic function, 
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with the fit parameters ci, ca, and C3. This monotonic func- 
tion describes the numerical data very well and vanishes for a 
completely detached particle, see Figs. [3] and [4] However, it 
cannot capture the partially-wrapped state encountered dur- 
ing the unwrapping transition at almost complete wrapping. 
This high wrapping fraction partially-wrapped state that has 
been found in Ref. is a very shallow energy minimum, 
which we cannot identify in our numerical calculations. Fur- 
ther analysis is done using the fit function, in particular the 
wrapping energy at any reduced adhesion strength is 
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Figs. [3] (a), § (a), and[4](c) show energies as function of the 
wrapping fraction at reduced tension B = 1 for wrapping of 
spherical, prolate, and oblate particles respectively. The open 
circles are the numerically calculated deformation energies 
that are fit by the solid line given by Eq. Q . The deformation 
energy thus increases monotonic with the wrapping fraction, 
the unwrapped state is stable. For finite adhesion strengths, 
the wrapping energies that are calculated using Eq. jHJ are 
non-monotonic functions of the wrapping fraction, such that 
partially and fully-wrapped states can be the stable states. 



III. 



RESULTS 



A. Wrapping energy 

The numerical data for the deformation energy without 
any adhesion in Figs. [3| (a), [4] (a), and [4] (c) is fit by the 
function given in Eq. ffl. For increased adhesion strength, 
the onset of wrapping occurs for the adhesion strength w\ 
for that the bending-energy cost equals the adhesion-energy 
gain at the contact point, see appendix [B] For all adhesion 
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FIG. 3. (a) Energies for wrapping a spherical particle as func- 
tion of the wrapping fraction A&d/A for reduced membrane 
tension 5 = 1. The figure shows the wrapping energy pro- 
files for six adhesion strengths: the numerically calculated 
data for zero adhesion strength and the corresponding fit 
function, E with equal energy for the non-wrapped and the 
completely- wrapped states (w = 5.99), the binding transi- 
tion Wi between the unwrapped and the partially-wrapped 
(w = 4.00), the binodal W2 between the partially- wrapped 
and the completely-wrapped state (w = 6.12), and the spin- 
odals S21 and S22 that are associated with Wi (for w — 2.46 
and w — 8.02 respectively). The phase boundaries separate 
5 regimes in the phase diagram with stable and metastable 
completely wrapped (CW), partially- wrapped (PW), and and 
non- wrapped states (NW); the stable state is underlined. The 
wrapping fractions that are plotted in Figj5] are marked by 
circles and the energy barrier shown in Fig. Mis indicated, (b) 
Sketches for spherical particles for wrapping fractions 0.25, 
0.50, 0.85 and 0.96, with the adhered membrane in red and 
part of the free membrane in blue. 



strengths w < wi the stable state for the particle is the un- 
wrapped state. Increasing the adhesion strength further, we 
find the envelopment transition from the partially-wrapped 
to the completely- wrapped state at adhesion strength u>2- 
For adhesion strengths wi < w < W2, a global minimum 
of the wrapping energy is found for a small wrapping frac- 
tion, < A a d/A -C 1. For adhesion strengths w > W2, a 
stable completely-wrapped state is found. The line E, where 
the energy of the unwrapped state equals the energy of the 
completely-wrapped state is a good estimate for the binodal 
W 2 . 

Whereas there is a continuous transition from the unbound 
to the bound state, the envelopment transition between the 
partially-wrapped and completely-wrapped state is character- 
ized by an energy barrier, AEbarrior- For W\ < w < w 2 , 
in addition to the stable partially-wrapped state at small 
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FIG. 4. Wrapping energies for (a) a prolate ellipsoid (PE) 
with aspect ratio 2 and (c) an oblate ellipsoid (OE) with as- 
pect ratio 0.75, analogous to Fig. [3] for a sphere, (b) and (d) 
Sketches for ellipsoids with wrapping fractions 0.30, 0.49, 0.84 
and 0.95. 



wrapping fraction, a metastable completely-wrapped state is 
found; similarly at higher adhesion strengths, in addition to 
the completely wrapped state a metastable partially-wrapped 
state is found. Indicated by the spinodal 5*22 that is associated 
with W2, the energy barrier between the metastable partially- 
wrapped state and the completely-wrapped state vanishes at 
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FIG. 5. Wrapping fractions for special points marked on the 
wrapping-energy functions in Figs. [3] and [4] for several val- 
ues of the reduced surface tension of the membrane and for 
several aspect ratios. At adhesion strength w 2 , the partially- 
wrapped state is found at wrapping fraction p 1 and the energy 
maximum at wrapping fraction p 2 ; P3 is the wrapping frac- 
tion for that the energy barrier between partially-wrapped 
and completely-wrapped state vanishes. 



an adhesion strength that is even larger than w 2 . Starting 
from a completely-wrapped state and continuously decreasing 
the adhesion strength, below a threshold value Wi a sponta- 
neous transition between the completely-wrapped state and 
the unwrapped state is observed, which is indicated by the sec- 
ond spinodal for W2, S21 [59]. The system shows strong hys- 
teresis, such that no partially wrapped state with small wrap- 
ping fraction is encountered during this unwrapping transi- 
tion. The height of the energy barrier between the partially- 
wrapped and the completely- wrapped state, ASbarrier, can 
be characterized by barrier height for w — W2, as indicated 
in Figs. [3] and [4] High energy barriers for ellipsoids with 
increased aspect ratio lead to an increased stability of the 
partially- wrapped state, compare section [III C| 

Special wrapping fractions are marked by points in Figs. [3] 
and|4j and are plotted in Fig. [5] as function of the surface ten- 
sion for several aspect ratios. The wrapping fraction pi de- 
notes the wrapping fraction for the stable partially-wrapped 
state, while P2 is the wrapping fraction for the energy maxi- 
mum, both for adhesion strength W2 ■ The saddle point of the 
energy, when the barrier between the partially-wrapped and 
the completely- wrapped state vanishes on the spinodal S22, 
occurs at wrapping fraction p$. All wrapping fractions for 
special points depend only weakly on the surface tension, but 
decrease slightly with the increasing surface tension at small 
surface tensions. 

When the partially-wrapped state has the same energy as 
the completely-wrapped state, the wrapping fractions for the 
partially-wrapped state are always well below 0.5. However, 
a strong dependence of the wrapping fraction on the particle 
shape is observed: whereas for prolate ellipsoids and spheres 
Pi « 0.2, for oblate ellipsoids it increases strongly with de- 
creasing aspect ratio. The maxima of the energy barriers are 
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FIG. 6. Wrapping states for spherical particles with reduced 
adhesion strength w and reduced surface tension a. Non- 
wrapped (NW), partially- wrapped (PW), and completely- 
wrapped (CW) states are found, stable and metastable states 
are indicated, stable states are underlined. The binodal 
Wi separates a stable non-wrapped and a stable partially- 
wrapped state, the binodal W% separates a stable partially- 
wrapped and a stable completely-wrapped state. The two 
spinodals S%\ and S22 belong to W%. For all states on E, 
the unwrapped state has the same energy as the completely- 
wrapped state. 



found for wrapping fractions between 0.6 and 0.7. For adhe- 
sion strengths w > W2 the partially-wrapped state becomes 
metastable; in this regime, the wrapping fraction at which the 
metastable partially-wrapped state appears increases with in- 
creasing adhesion strength until it reaches the critical value 
P3, which is found at approximately half wrapping. 



B. Wrapping diagrams 
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FIG. 7. Wrapping states for prolate ellipsoidal particles with 
aspect ratios (a) 1.5 and (b) 2, plotted analogously to Fig. [6] 



For wrapping a spherical particle by a tension-less mem- 
brane, the only relevant energies that determine the uptake 
process are bending and adhesion energy for the adhered part. 
At every point of the sphere the bending energy is constant 
and the free membrane around the particle forms a catenoid- 
like deformation, i. e., a minimal surface with vanishing bend- 
ing energy cost. Thus there exist only 2 possible stable states, 
viz. a non- wrapped and a completely-wrapped state, that are 
separated by a continuous transition at w — 4. For a spherical 
particle and a finite membrane tension, also a stable partially- 
wrapped state is found, compare section [ill A| For ellipsoidal 
particles this partially- wrapped state exists both for finite and 
for vanishing surface tension. 

In the wrapping diagrams for surface tension and adhesion 
strength in Figs.[3J[7] and[3] five regimes can be identified with 
different combinations of stable and metastable unwrapped, 
partially wrapped, and completely wrapped states. For small 
adhesion strengths, a stable non-wrapped state is found. In 
between the spinodal for the spontaneous transition between 
the completely- wrapped and the non- wrapped state S21, and 



the binodal for the binding of the colloid to the membrane 
Wi , in addition to the stable non- wrapped state a metastable 
completely-wrapped state with high wrapping fraction ap- 
pears [50] . The binding transition occurs at adhesion strength 

w>i,pe = [l + (6/a)~ 2 ] and mj^oe = 4(6/a) 2 for prolate and 
oblate ellipsoids respectively and is independent of the surface 
tension, see appendix |B| and Refs. 30, 35, and 61. Beyond the 
threshold adhesion strength for binding, a partially-wrapped 
state coexists with the metastable completely-wrapped state. 
For adhesion strengths beyond those for the binodal W2, the 
completely-wrapped state becomes stable and the partially- 
wrapped state becomes metastable. Finally, beyond the spin- 
odal 522, the energy barrier between the completely- wrapped 
and the partially-wrapped state vanishes. 

For spherical particles, all phase boundaries meet at a triple 
point for zero surface tension, see Fig. [6] However, as shown 
in Figs. [7] and [8] for ellipsoids, even for zero surface tension 
Wi is shifted to smaller values of the adhesion strength com- 
pared with spherical particles, while the binodal W% is shifted 
to higher values of the adhesion strength. For comparable as- 
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FIG. 8. Wrapping states for oblate ellipsoidal particles with 
aspect ratios (a) 0.75 and (b) 0.5, plotted analogously to 
Fig.© 
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FIG. 9. Energy barriers between partially-wrapped and com- 
pletely wrapped state for the adhesion strength W2 where par- 
tially and completely-wrapped state have equal energies. The 
barrier energies are plotted for several aspect ratios as func- 
tion of the reduced surface tension a and fit by a function 
A-Bbarricr = Eq + E\b v . The values for the sphere are com- 
pared with those from Ref. 1351 
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pect ratios, this effect is stronger for oblate ellipsoids than for 
prolate ellipsoids. This range of adhesion strengths, for which 
stable partially-wrapped states are found, increases both with 
the aspect ratio of the particle and the surface tension of the 
membrane. For increased aspect ratio it is easier to attach 
particles to the membrane, but at the same time it is more 
difficult to achieve completely-wrapped state. Binding oc- 
curs already for smaller adhesion strengths because of the re- 
duced particle curvature at the point of first contact, whereas 
complete wrapping requires higher energies because of the in- 
creased curvature of the particle at the tips. For high aspect 
ratios and rather 'flat' particles, the dependence of the bin- 
odal W2 and both associated spinodals on the surface tension 
decreases. 



FIG. 10. Energy barrier Eo as a function of the aspect ratio 
b/a for vanishing surface tension. Oblate particles have b/a 
< 1 and prolate particles have b/a > 1. 



C. Energy barrier 

In Figs.[3]and[4] the wrapping energy for adhesion strength 
u>2 shows a partially-wrapped state at a wrapping fraction 
pi that is separated by an energy barrier Awarder from the 
completely- wrapped state. For w > W2, the height of the 
barrier between a metastable partially- wrapped state and the 
completely-wrapped state determines the dynamics for the 
wrapping process. The activation energy that is needed for 
complete wrapping can be provided either by thermal kicks 
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or by active forces from motor proteins. A comparison of the 
height of the energy barrier shown in Fig. [9] with the thermal 
energy ksT therefore allows to estimate a characteristic time 
for the transition |13j . 

The energy barriers vary with the surface tension of the 
membrane and with particle shape. While for spherical par- 
ticles the energy barrier vanishes for the tension-free case and 
can be fit by a power law [35], it levels off to finite values both 
for prolate and oblate ellipsoids. The energy barriers for ellip- 
soidal particles thus can be orders of magnitude higher than 
those for spherical particles. We fit the energy barriers as 
function of the wrapping fraction by AZ?barricr = Eo + E\a v . 
While the exponent v is similar for all particles, the prefactors 
Eo and E\ strongly depend on the particle shape. 

While in section |IH B| an extended parameter range for a 
partially-wrapped regime has been found for ellipsoidal par- 
ticles compared with spherical particles, the energy barrier 
characterizes the stability of metastable partially-wrapped 
states. This can practically imply that partially-wrapped el- 
lipsoidal particles are found in experiments even for adhesion 
strengths much higher than u)2, for that wrapping diagrams 
already predict a completely-wrapped particle. In Fig. |I0| 
the energy barriers for vanishing surface tension are plot- 
ted as function of the aspect ratio of the particle. Only for 
Eo < k&TKnK) the barrier height is comparable to thermal 
energy and complete wrapping is expected to occur at 1W2. 
Eo is of the order of 1 already for prolate ellipsoids with as- 
pect ratio 1.5 thus corresponds to an energy barrier of about 
7tk. This implies that for particles whose shape considerably 
deviates from a sphere, adhesion strengths that are consid- 
erably higher than wi are required to wrap the particle and 
the partially-wrapped regime might extend almost up to the 
spinodal Sii- 



D. Role of shape anisotropy 

Stable non-wrapped, partially-wrapped and completely- 
wrapped regimes can be identified in Fig. 1 1 1 1 where the state 
of the system depends on membrane surface tension, parti- 
cle aspect ratio, and adhesion strength. For spherical parti- 
cles, the transition from the non-wrapped to the completely- 
wrapped state occurs at the triple point for w — 4. For ellip- 
soids, the transition from the non-wrapped to the partially- 
wrapped state is found already for Wi < 4, where ibije — > 1 
for prolate ellipsoids with high aspect ratio and IUi,oe — > 
for disc-like oblate ellipsoids, compare Eqs. (B3l and (B4l. 



The limiting cases for wi correspond to an infinite cylinder 
[371 138] and a flat disc [451 162] respectively. For the envelop- 
ment transition, we find an increased adhesion strength ibi 
for the prolate and oblate ellipsoids with increasing aspheric- 
ity; an additional shift towards higher adhesion strength is 
observed for finite surface tension. 

Besides the numerical minimization technique to calculate 
W2, we have used a deformed-catenoid approximation to esti- 
mate the wrapping energy that allows to estimate W2 even for 
high aspect ratios such as 10. For a tensionless membrane, we 
numerically calculate the deformation energy for a wrapped 
ellipsoid by deforming a sphere with a catenoidal membrane 
patch, without any actual minimization. This method works 
well, because without surface tension the contribution of the 
deformation of the membrane patch around the particle is 
small. For oblate particles and bending-only, the membrane 
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FIG. 11. Stable wrapping states for different aspect ratios 
b/a, reduced adhesion strengths w, and for a = and for 
a — 2. Prolate ellipsoids (PE) have aspect ratios b/a > 1, 
while for oblate ellipsoids (OE) have aspect ratios b/a < 1. 
Stable wrapping states are the non- wrapped state (NW), the 
partially- wrapped state (PW), and the completely- wrapped 
state (CW); the roman numbers correspond to those used in 
Figs. [6] [7] and [8] For blue open squares labeled with index 
dc, the deformed catenoid approximation has been used. 



patch around the particle will assume a catenoid shape with 
vanishing bending energy and therefore does not have to be 
calculated. 

As discussed in section |III C[ for w >Wi the stable state 
is the completely- wrapped state ( CW ), but still the particles 
may be stuck in metastable partially-wrapped states unless 
there are sufficiently large fluctuations such that the energy 
barrier can be crossed, compare Fig. [9] Thus, for a given ad- 
hesion strength that is smaller than the one for the spinodal 
for spontaneous wrapping, we predict a considerably smaller 
amount of completely-wrapped particles the stronger the as- 
phericity is. 



IV. SUMMARY & CONCLUSIONS 

Complete wrapping a spherical particle using a tensionless 
free membrane occurs at reduced adhesion strength w — 4, di- 
rectly from the unwrapped to the wrapped state without any 
energy barrier. For a membrane with surface tension the en- 
velopment transition is shifted to higher adhesion strengths, 
while the binding transition remains at w — 4. The new 
partially-wrapped state is separated from the completely- 
wrapped state by an energy barrier. We find that a partially- 
wrapped state also exists for ellipsoidal particles, in a wider 
region of the phase diagram and with a higher energy bar- 
rier to the fully-wrapped state than for spherical particles. 
Therefore, the spherical shape facilitates complete wrapping 
for single particles; a biological example is the uptake and 
budding of almost spherical viruses [631 164] . 

In addition to a shift due to a membrane surface tension, 
for ellipsoidal particles the envelopment transition addition- 



ally shifts to higher adhesion strengths with increasing non- 
sphericity. The binding transition shifts to smaller adhesion 
strengths: for very long prolate ellipsoids w% — ¥ 1 and for 
very flat oblate ellipsoids to wi —¥ 0. The partially-wrapped 
state is additionally stabilized by a higher energy barrier to 
the completely-wrapped state. Therefore, attachment of el- 
lipsoidal particles to a membrane and partial wrapping is 
facilitated compared with spheres, while complete wrapping 
is hindered; typical wrapping fractions are 20 — 40%. Elon- 
gated viruses are found to form patterns on the membrane, 
see Rcf. 25 However, similar to curved inclusions the parti- 
cles may bud cooperatively [91 112j. 

For typical lipid bilayer bending rigidities of k = 20 ksT, 
the reduced energies E ~ 1 correspond to typical energies 
E ~ 60 ksT. Energy barriers at adhesion strengths where 
the energy of the partially-wrapped state equals the energy 
of the completely wrapped state can thus be of the order of 
100 ksT. Often the length scale for the particle in soft matter 
and biological system is of the order of hundred nanometer, 
therefore a = (a 2 / ' k)o = 500 a nm 2 /ksT; a biologically rele- 
vant surface tension for the cell membrane, a = 0.003 dyn/cm 
|65| corresponds to a « 1 for a particle size a = 100 nm, to 
a w 0.25 for a particle size a — 50 nm, and to a w 4 for a 
particle size a — 200 nm. At a reduced tension a — 0.25, the 
energy barrier to the completely wrapped state is an order 
of magnitude higher for prolate ellipsoids of aspect ratio 2 
and oblate ellipsoids of aspect ratio 0.5 than for a spherical 
particle. 

A typical adhesion strength can be estimated based on the 
binding strength of the HIV virus, w « O.lfcBT/nm 2 [63] . 
The reduced adhesion strength is thus w = 2(a 2 / 'k)w = 100 
for a particle with size a — 100 nm, w ~ 25 for a — 50 nm, and 
to w rs 400 for a = 200 nm. All adhesion strengths are well 
in the region where wrapping occurs. However smaller adhe- 
sion strengths may occur for other viruses or smaller recep- 
tor densities in the membrane. An increased reduced tension 
hinders budding, but because a and w both scale with the 
squared size of the particle, according to our calculations and 
previous wrapping calculations larger nano-particles are more 
likely to be wrapped completely. A lower limit for uptake is 
« 20 nm, when the adhesion energy balances the bending en- 
ergy 3, 66 68 . However, there is also an upper limit given 
either by the length scale \fnfo where the wrapping becomes 
surface-tension dominated [35] or by receptor availability |66f - 
[68]. 

Our theoretical calculations predict an enhanced stability 
of partially-wrapped states for ellipsoidal particles. For cells, 
a lower uptake for ellipsoidal particles has been found experi- 
mentally [161 17 20 . partially adhered ellipsoidal particles are 
also discussed in particular in Ref. 1201 similarly disk-shaped 
particles show an increased adhesion to membrane and lower 
uptake [691 170) . Receptor-mediated wrapping of ellipsoidal 
particles has been studied in Ref. [71] with similar findings as in 
our work, but without surface tension of the membrane. The 
partially wrapped state can be of advantageous both from an 
application point of view as well as from a biological point of 
view: for example, elongated particles can be used as markers 
for imaging that stay in the cell membrane 27, 28 and Ebola 
and Marburg viruses are not easily taken up by macrophages 
and may thus have a high virulence [3]. 
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Appendix A: Triangulated membranes 

Triangulated membranes are a powerful tool to study mem- 
branes and interfaces 5,8,49 54. Using Surface Evolver [49]. 
the membrane shape can be minimized with different schemes 
and the triangulation can be refined at any stage. In this ap- 
pendix, we present the discretization that has been used to 
calculate the energy in Eq. |2|. Surface tension and adhe- 
sion energies are proportional to the area of the membrane, 
they can be calculated basically as sum over all triangle areas. 
The bending energy is calculated, using the Surface Evolver 
method "star_perp_sq_mean_curvature" , which assumes every 
vertex has incident triangles forming a star network around 
it, such that a v is the average area associated locally with the 
vertex. The force acting on a vertex when the area changes 
thereby causing the vertex to move is given by the gradient 
of the area (Va„) associated with the vertex. For a smooth 
surface the gradient of the volume (VV V ) is equal to the mag- 
nitude of the area. But for a triangulated patch, the facets 
around the vertex tilt and thus the area is greater than the 
magnitude of the gradient of the volume. For a membrane 
patch, with a central vertex v and neighboring vertices v 1 , 
v 2 , v 3 , ... v n a volume is given by 

V v = -v ■ [v 1 x v 2 + v 2 x v 3 + ...v n x v 1 ] , (Al) 
6 

and the gradient of the volume is therefore 

VV„ 4k'xv ! + v a xv 3 + ...v n x v 1 ] . (A2) 
6 

The local mean curvature (h v ) at the vertex is 

fe " = 2 VV„-VV„ - (A3) 

Thus the discretized form of total squared mean curvature 
integral is given by, 

n 

£bend = 2ft ^2 a vhl ■ (A4) 

U = l 

Once assembled, the surface can be minimized using 
different schemes, default being moving towards the direction 
of steepest descent of the energy linearly, while the mesh 
may be refined or smoothed at any stage. One may employ 
other minimization schemes like Hessian approach which 
calculates the energy of the surface configuration for a 
small perturbation and then uses the Hessian, a square 
matrix formed of the second derivatives of the energy which 
determines the best quadratic approximation of the energy 
to look for minimum energy states. Once converged to a 
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minimum energy state, the surface may be analyzed for 
deformation profiles. 



Appendix B: Binding transition (Wi) 

The binding transition for particles to a membrane is de- 
termined by the mean curvature of the particle at the contact 
point and is independent of surface tension 35 . To calcu- 
late the binding transition for any particle shape, the Monge 
parametrization can be used, where the surface is described 
by a height field, h(p), where p = (x,y) are the coordinates 
in the reference plane. For an almost planar membrane, a 
small-gradient approximation for the bending energy gives 



£def = J dA^{V 2 h{ P )f , 
with J dA the integral over the reference plane. 



(Bl) 



The critical adhesion strength at the binding transition 
marks the onset of adhesion, separating the unwrapped 
regime from the (partially) wrapped regime. Because it is 
completely determined by the competition between adhesion 
energy and bending energy, at the transition the mean curva- 
ture must equal the adhesion strength wi, 



wi = a 



[v 2 Mp = o)] 2 



(B2) 



for the contact point between membrane and particle at p = 
0. For spheres [\7 2 h(p — 0)] 2 = 4/a 2 , we find w 1 — 4, for 
prolate ellipsoids of aspect ratio b/a, such that the they are 
wrapped with their long axis parallel to the membrane, 



Wi(PE) = 



and for oblate ellipsoids, 



1 + 



wi(OE) = 4(6/o) 



(B3) 



(B4) 
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